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Abstract 

We investigate the observables available in the angular distribution of -B — t- K* to identify 

those suitable for measurements in the first few years of LHC data taking. As experimental 
uncertainties will dominate, we focus on observables that are simple to measure, while maxim- 
izing the potential for discovery. There are three observables that may be extracted by counting 
signal events as a function of one or two decay angles and correspond to large features of the full 
angular distribution in the Standard Model: A-p-Q, Fl, and S5. Two of these are well known in 
the experimental community; however, we show that measuring 5*5 adds complementary sens- 
itivity to physics beyond the Standard model. Like ^fBj it features a zero-crossing point with 
reduced hadronic uncertainties at leading order and in the large recoil limit. We explore the 
experimental sensitivity to this point at LHC 6 and show that it may be measured with high 
precision due to the steepness of the distribution. Current experimental model independent 
constraints on parameter space are presented and predictions made for the values of the Ape 
and 5*5 zero-crossing points. The relative impact of LHC 6 measurements of Ayb-, Fi^, and 55, 
with 2fb~^ of integrated luminosity, is assessed. These issues are explored with a new model of 
the decay that can be used with standard simulation tools such as EvtGen. 



Keywords : S-Physics; Beyond Standard Model; Rare Decays. 



* a.k.m.bharucha@durham.ac.uk 
twill. reece@cern.ch 



1 Introduction 



The decay — t- K*"fi~^fj,~ is a golden channel for the study of flavour changing neutral currents 
(FCNC) at the Large Hadron Collider (LHC). The four-body final state, as K*^ — )• Ktt, means 
that there is a wealth of information in the full-angular distribution that is complementary 
to that available in the widely studied 6 — )• 57 decays. In the presence of physics beyond 
the Standard Model (SM), new heavy degrees of freedom may enter the 6 — )• s loops. These 
can alter the decay amplitudes, affecting the full-angular distribution observed. This makes 
Bii — )■ K*^fi^fi^ one of the most promising places in the flavour sector to search for new physics 
(NP) at the LHC (see Ref. [1] for a review). We concentrate on the large-recoil regime, where 
the energy of the K* is large such that QCD factorization is applicable. The low-recoil regime 
was described in Ref. |2], however at present form factors in this regime are not well known. 
A number of interesting measurements have already been made [3HS]. They are broadly in 
agreement with SM predictions; however, experimental precision is currently too low for firm 
conclusions to be drawn. 

The properties of the full-angular distribution have been studied by many authors and a 
number of potential measurements have been identified; e.g. Refs |lUH16j . Particular emphasis 
has been placed on finding angular observables with reduced theoretical uncertainties or en- 
hanced sensitivity to particular classes of NP. However, in the first few years of LHC data 
taking the dominant sources of uncertainty will be experimental; thus, the emphasis should be 
on finding quantities that can be cleanly measured with relatively small uncertainties. Once 
very large data sets have been collected, it will be possible to use a full-angular analysis to 
extract the various underlying amplitudes directly \13\ [T7] . This will allow the determination 
of many theoretically clean observables. However, performing this kind of analysis will not be 
possible until detectors are very well understood and the number of collected signal events are 
in the thousands. Prior to this, symmetries and asymmetries of the full-angular distribution 
can be used to extract some observables individually from angular projections |14 ^ 1151 [T8ll20j . 

In this paper, we focus on observables that correspond to large features in the — t- 
K*^ IJL^ full-angular distribution and can be measured by counting the number of signal events 
as a function of one or two decay angles. We then investigate the relative experimental sensit- 
ivities to these observables at LHC 6 [21] and their projected impact on the allowed parameter 
space after measurements with 2 fb~^ of integrated luminosity. The rest of the paper is struc- 
tured as follows: In the next section we give a brief overview of the theoretical framework 
employed with details of the decay amplitude calculation; in Sec. [3j observables that will be rel- 
evant for analyses with the first few years of LHC data are discussed, and details of benchmark 
NP models provided. We also summarize the impact of existing experimental measurements on 
constraining the NP contribution to the Wilson coefficients. In Sec. [5| we analyse the possibility 
of detecting NP effects at LHC 6 using our chosen observables. In Sec. |6| the potential impact 
of these measurements on parameter space is assessed. Finally, in Sec. [7j a short summary is 
given. 
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Ci(m) 






C4(m) 






-0.135 


1.054 


0.012 


-0.033 


0.009 


-0.039 






ACf (^) 






-0.306 


-0.159 


4.220 


-4.093 



Table 1: SM Wilson coefficients at = m,, = 4.52GeV/c2, where /^Cf{^x) ~ Cf{^x) - Y{q^). 



2 Theoretical Details 
2.1 Introduction 

A decay model following Ref. |22j has become the standard tool for studies of — )• K*^fi'^n~ 
within the experimental community due to its inclusion in the decay simulator EvtGen [23]. A 
significantly improved version of that model with much greater support for the simulation of NP 
as well as a state-of-the-art SM treatment has been developed as part of the present work [24] . 
We present our theoretical framework in a way that allows direct comparison with Ref. [22], 
by expressing the decay amplitude in terms of the auxiliary functions used in that reference. 
Calculation of these requires Wilson coefficients, form factors and quantum-chromodynamics 
factorization (QCDF) corrections, as described in detail in this section. 



2.2 Wilson Coefficients 

The Wilson coefficients, Ci{fi), are process-independent coupling constants for the basis of 
effective vertices described by local operators, Oi{^), and encode contributions at scales above 
the renormalization scale, fi. For a given NP model, new diagrams will become relevant and the 
Cj(//)'s may change from their SM values; additional operators may also become importanllM 

(u) 

The weak effective Hamiltonian, neglecting doubly Cabibbo-suppressed contributions, "Hg^ , is 
given by 

^eff = -^Ai |^CiOj + C20^ + J^ao, + ^(c,o, + c;.o;.) j , (i) 

where j = 7, 8, 9, 10, P, S, Gp is the Fermi constant, and At = VtbV^* is the relevant combination 
of Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. The operators O and O' are defined 
in Ref. [15], and a subset is given explicitly in App. [Aj 

The primed operators have opposite chirality to the unprimed ones and their correspond- 
ing coefficients, C-(/i), are suppressed by nis/mh or vanish in the SM; however, they may be 
enhanced by NP. We neglect the contributions from O'^ for 1 < i < 6 as they are either heavily 
constrained by experimental results or generically small; NP contributions to Oj_iq may still 
be important and are included. We also include the scalar and pseudoscalar operators Cgp- 
These vanish in the SM but may arise in certain NP scenarios, for example in the case of an 
additional Higgs doublet. 

The Wilson coefficients are calculated by matching the full and effective theories at the scale 
of the W boson mass, m^y. For the SM Wilson coefficients, we aim at next-to-next-to- leading 



comprehensive review of efTective field tlieories in weak decays can be found in Ref. |25) 
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logarithmic (NNLL) accuracy. This requires calculating the matching conditions at /i = mw to 
two-loop accuracy. This has been done in Ref. [26]. NP contributions are included to one-loop 
accuracy only, as two-loop corrections are expected to be small. This was shown explicitly for 
the MSSM in Ref. j27j . The Wilson coefficients must then be evolved down to the scale fi ~ ruh. 
The evolution has been implemented using the full 10 x 10 anomalous dimension matrix following 
Refs [28ti30] . The primed operators, Oj_iq, are evolved as their unprimed equivalents; however, 
the scalar and pseudoscalar operators O^jp are defined to be conserved currents and do not 
mix with the other operators and so do not require evolution. For convenience, we define the 
following combinations of Wilson coefficients: 



in 1 4 20 80 
a= 3 9 3 9 



Cf^ — ^8 + C3 — - C4 + 2OC5 — Cq] 

tts 3 

Cf = -C + Yiq^); 
as 

r<cS _ . 
'^10 — — ^10) 
as 



f,/eS - ■ (9) 

•-"7,8,9,10 — ^ "-^r.s.g.io 5 v^^ 

Us 



where is the invariant mass squared on the muon pair and Y{q^) is defined in Ref. [Slj- Tab.|l| 
gives the values of the Wilson coefficients at /i = mb_ps(2GeV) in the SM. The treatment of 
quark masses in the PS scheme is discussed in Sec. |2.6[ 



2.3 Form Factors 

Bd — )• K*^fi~^fi~ is characterized by eight form factors, V{q'^), ^0-3(9^) and Ti_3(g^). These 
are hadronic quantities that, for certain ranges in g^, may be obtained by non-perturbative 
methods. Their definition in terms of hadronic matrix elements can be found, for example, in 
Ref. [3^. Lattice field theory currently offers a prediction for the form factor Ti(0) relevant to 
B — )■ K*'y [33], but not for the others. However, QCD sum rules on the light cone (LCSR) is a 
well established alternative technique that provides results for the desired range in [151 [32] . It 
is an extension of classic QCD sum rules [34j, in which matrix elements are evaluated via both 
operator product expansion and dispersive representation. Quark-hadron duality then leads 
to sum rules for the desired hadronic quantities. LCSR follows a similar procedure to obtain 
sum rules for the form factors, but the operator product expansion in terms of vacuum conden- 
states is replaced by a light-cone expansion in terms of universal light-cone meson distribution 
amplitudes. A comprehensive review of QCD sum rules and LSCR can be found in Ref. |35| . 

We use the full set of LCSR form factors in our model \32\ 136] , where the sum rules for all 
form factors except for ^0 were calculated at 0{as) accuracy for twist-2 and-3 and tree-level 
accuracy for twist-4 contributions. Note that the normalization of the form factors we use differs 
slightly from Ref. |15| . however this will not have much impact on the observables, as they are 
normalized by the total decay rate, so the effect will cancel out. We estimate the uncertainties 
using the values provided in Ref. [32] for = 0, as shown in Tab. I2j Note that ^3(0) and 
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m 


Atot 




V 


0.411 


0.033 


0.44 5a, 


Ao 


0.374 


0.034 


0.39 6a, 


Ai 


0.292 


0.028 


0.33 6a, 


A2 


0.259 


0.027 


0.31 6a, 


Ti 


0.333 


0.028 


0.34 6a, 


T3 


0.202 


0.018 


0.18 6a, 



2.3 



Table 2: Form factors for — ?• K iT' from LCSR at g = [32j, as described in Sec. 
Here Atot is the total error arising from the uncertainty on all input parameters with the 
exception of the Gegenbauer moment ai. A^^ contains the uncertainty due to ai, where 6a, is 
defined 6a, = ai{K*, IGeV) - 0.1. 



T2(0) are not included in the table, as they can be found using the relations ^43(0) = Aq{0) and 

r2(o) = ri(o). 

In the large energy limit of the K* , the form factors satisfy certain relations and, therefore, 
can be reduced to two heavy-to-light or soft form factors, denoted ^± and [STHIU] . These 
reduced form factors are generally used within the QCDF framework |31|, IHj . The relations 
are studied through appropriate ratios of the LCSR predictions for the full form factors in 
Appendix B of Ref. [15]. It is shown that those involving ^± are almost independent of q'^, but 
those involving ,^|| have a definite dependence on q"^, so are probably more sensitive to the l/rrih 
corrections neglected in QCDF. 



2.4 QCD Factorization Corrections 

QCD factorization is a framework in which the 0{as) corrections to Bd — t- K*^ijl^ ii~ can be 
calculated in the combined heavy-quark and large-recoil energy limit; this applies when the 
energy of the K* is large. These corrections take into account contributions that cannot be 
included in the form factors, such as the non-factorizable scattering effects arising from hard 
gluon exchange between the constituents of the B meson. 

Our calculation of the decay amplitude includes QCDF corrections at next-to-leading order 
(NLO) in as but leading order (LO) in l/mf,. These corrections are included in the definitions 
of T\\{q^) and T±{q'^) found in Ref. [3l] and are given in terms of and ^y; however, 0{as) 
factorizable corrections that arise from expressing the full form factors in terms of and 
must then be subsumed. Following Ref. [I5], we instead express our LO results for the decay 
amplitude in terms of the full form factors. Factorizable corrections are then redundant and 
the main source of 0{l/mb) corrections is automatically included. In addition, we neglect weak 
annihilation corrections at LO in 1/mf, and 0{as) as they are dependent on the numerically 
small Wilson coefficients C3 and C4. 

We denote 7jj^^°(g^) and Tf^'^iq^) to be the analogues of r\\{q^) and 71 (g^) from Ref. [3l] 
with the only relevant 0{as) contributions included. We also define 7jj''^'"^((j'^) and Tf^^'~'{q'^); 
the primes indicate that the unprimed Wilson coefficients should be replaced by their primed 
equivalents. In order to extend the results of Ref. to include NLO corrections, we must 
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make the following replacements: 

TUB 

where EK*{q^) is the energy of the K* and is the mass of the B meson. 

We have now introduced the Wilson coefficients, form factors and defined the QCD factor- 
ization corrections. These are all ingredients for the auxiliary functions describing the decay 
amplitude, as seen in the following subsection. 



2.5 Decay Amplitude 

The Hamiltonian defined in Eq. ([T]), combined with the standard definitions of the form factors, 
leads to the following decay amplitude [22l 



[7;,i(/Z7V)+r^'(/"7^5^)+5(/x/i)] (4) 

where 

= A{q^)e^p^pe'PrBpK* " iB{q^) + iC{q^W ■ pb)p^ + iD{q''){e* ■ ps) % (5) 

= E{q^)e^p^p e*Pp%p^^, - iF{q^) e*^ + iG{q^)ie* ■ pb)p^ + iH{q^){e* ■ pb) q^ (6) 

and 

S = i2mK^{e* ■pB)I{q^). (7) 

Here, pb,k* and tub^k* are the four- momenta and masses of the respective particles in the 
B meson rest frame, p = pb +Pk*^ Q = Pb —PK*, and e* is the K* polarization vector. The 
circumflex denotes division by tub (e.g. itik* = ttik* /ms)- The auxiliary functions A-I{q^) 
follow Ref. [22]; however, we have updated the previous expressions to include additional primed, 
scalar, and pseudoscalar operators, as well as QCDF correction via ^^^^ {q^) and t[_ ^^^'^ {q^) 



as outlined in Sec. 2.4 They are defined as: 



A{q^) =--^{Cf + C',^^)Viq') + ^ ( (Cf + C'.^^mq'^ 
1 -I- rriK* q 

B{q') =(1 + mK^)l^{Cf - Cr)^i(^') + ^(1 -rhK'){ {Cf - C'.^^mq'' 

+ 2EK^{q^){T^^^{q^) - T^^Olg^)) ) ); (8b) 
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Parameter 


T 7" 1 

Value 


T~» r 

Rer. 


rus 


0.104 


[43J 


mc,Ps(0.7GeV) 


1.5 GeV 


[44J 


m6,ps(2GeV) 


4.52 GeV 


[45J 


rhtirht) 


162.3 GeV 


146J 



Table 3: Quark masses 



1 - rhj^^ 

+ 2m,((C7f - C'.^^miq') + "^^Uq')) 
+ (1 + (l-^l^*y£^^*(g') )(7-NLO(^2) _ 7-;NL0(^2)) 
+ 7^NLO(^2)_^/NLO(^2) 

^(«') =7TT^(^io + Clf (8d) 

(i + rriK* ) 

F{q') =(1 + m^.)(Qt - Clf (8e) 
G(g') ={Cll - C[f)-^f^y (8f) 



^2(g2) 



1 



H{q') =^(Cit - ) ( (1 + rhKMiiq') - (1 - ^K*)A2{q') 



- 2mK^-A,{q')] - '^^f^A,{q^){Cp - C'p); (8g) 
I{q^)=-Ao{q^){Cs-C's). (8h) 
The recoil energy of the K* is given by 

2x _ m| + m\, - q^ 

• (9) 

Using the equations of motion for the muons, 

(fimtil^) = and q^{Ulfil5lA = -2m^/i75 /i, (10) 

where is the muon mass, we see that D{q^) vanishes and H(q^) is suppressed by a power of 
m^. However, H[q^) receives a pseudoscalar contribution inversely proportional to allowing 
for some sensitivity to Cp — C'p [3^. The observables described in Sec. 3.1 (e.g. Eqs (16)-(17)) 
may be calculated directly from the amplitudes given in Eq. Q; the necessary formulae are 
presented in App. [B] and implemented in our model. 
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Parameter 


Value 


Ref. 


Parameter 


Value 


Ref. 


fB 


200 ± 25 MeV 


m 


af j^, (2 GeV) 


0.03 ±0.03 


m 


Xb{2.2 GeV) 


0.51 ±0.12 GeV 


m 


al ^, (2 GeV) 


0.02 ±0.02 


m 


/i*(2 GeV) 


163 ± 8 MeV 


m 


ai^,{2 GeV) 


0.08 ±0.06 


m 


fl 


220 ± 5 MeV 


m 


(2 GeV) 


0.08 ±0.06 


m 



Table 4: Hadronic parameters 



2.6 Numerical Input 
2.6.1 Quark Masses 

The calculation of the auxiliary functions requires the bottom quark pole mass, which is known 
to contain large long-distance corrections. To avoid this, a renormalization scheme, known as 
the potential subtraction scheme (PS), was introduced in Ref. H?]. The quark mass defined in 
the PS scheme has the advantage that the large infrared contributions are absent, while being 
numerically close to the pole mass. It is suitable for calculations in which the quark is nearly 
on-shell. Following Ref. [3l], we replace the pole mass by the PS mass, mps(/i/), using 

m = mps(/i/) ± ^/i/ ± ©(a^) (11) 

and neglect any resulting terms of ©(a^). Here /i/ is the scale at which the PS mass is calculated. 
All occurrences of the symbol nn, in our formulae refer to the PS mass, mb^ps(2 GeV), as shown 
in Tab. [31 

The operator O'j is defined in terms of the modified minimal subtraction (MS) mass. In the 
MS scheme, the 1/e poles are simply removed, along with the associated terms in 7 and 47r. 
Therefore, when the h quark mass arises in combination with Cj^, we replace the MS mass, rh, 
by the pole mass, using 
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7fi(/i) = m (^1 ± ^ (^3 In -| - 4j ± 0(a2) J . (12) 
This leads to factorizable OioLg) corrections to T^YPif^^ and T'^pf^icj^) as found in Ref. [31] . 



For consistency, we calculate the charm quark pole mass using Eq. (11). Here the PS mass 
is taken from the most recent calculation as in Tab. [3| The resulting pole mass agrees with 
results in Ref. [l3], where it is calculated from the MS mass. The top quark mass enters the 



calculation of the Wilson coefficients, and for this we use the MS mass in Tab. 3l as in Ref. [T5i 



2.6.2 Hadronic Parameters 



In addition to the form factors described in Sec. 2.3, the QCDF corrections require light-cone 
distribution amplitudes and decay constants. The light-cone distribution amplitude for both 
the B and K* mesons enter the hard scattering corrections. For the B meson we follow the 
prescription in Ref. [31] using the values for Ab given in Tab. |4| For the K* meson we use the 
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T~» J. 

Parameter 


Value 


Parameter 


Value 


niB 
Mw 


5.28 GeV 
0.896 GeV 
0.106 GeV 
80.4 GeV 


Kb 

7 


0.226 ±0.002 
(3.93 ±0.36)10-3 
(77^30)0 
(1.166)10-5 GeV-2 



Table 5: CKM matrix parameters, additional masses and constants from Ref. |43| . 



standard Gegenabauer expansion, 

= 6n(l - n)(l ± a^^^.cf /2)(2'u - 1) ± a'^'j^,ci^^^\2u - 1)), (13) 

for m =±, II, taking the coefficients from Tab. |4j We also require the decay constants for both 
the B and K* mesons. Additional parameters are summarized in Tab. [5} 



3 Observables and New physics 

Having established the basic theoretical framework, we proceed to discuss experimental observ- 
ables for Bd K*V^/""- 



3.1 Observables 

The full-angular decay distribution can be written as: 

— = ^i{q^,ei,eK*,(p), (14) 

dq^ d cos di d cos c^k* dcp sZtt 

where the angles 9k* , Oi and (p are defined as follows: 9k* is the angle between the and B 
in the rest frame of the K* , and is defined in the range — 1 < cos 9k* < 1; 9i is defined as the 
angle between the //- and B in the di-muon centre of mass frame, and is defined in the range 
— 1 < cos 9 1 < 1; is the angle between the normal to the K-ir plane and the normal to the 
di-muon plane, and is defined in the range < i?^ < 27r. For the conjugate decay, the angles are 
defined analogously, but with reference to the and /i"*". We can then express /(g^, 9i, 9k*, (p) 
in terms of these angles as follows: 

/(g^ 9i, 9k* , </>) = If sin^ 9k* + If cos^ 9k* + (/| sin^ 9k* ± /| cos^ 9k* ) cos 29i 
± I3 sin^ 9k* sin^ 9i cos 2<f> + I4 sin 29k* sin 29i cos (p 
+ I5 sin 29k* sin 9i cos 4> 

+ (/g sin^ 9k* + Iq cos^ 9k* ) cos 9i + I7 sin 29k* sin 9i sin (p 

+ Ig sin 29k* sin 29i sin (p + Ig sin^ 9k* sin^ 9i sin 2(p. (15) 

The angular coefficients where z = 1 to 9 and a = s or c, describe the decay distribution. 
A natural set of observables was identified in Ref. [15] by taking combinations of these /j-"^'s 
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that emphasize CP-conserving and CP-violating effects. These were defined as 



/ffi^, (16) 

where the Al s have also been studied in Ref. [H]. We introduce the rate average, which, 
for a variable V{q'^), is given by 

iGeV^ 



(^>l-6GeV^= / , dq^V{q 
JlGeV^ V 




dq' 



2 • (18) 



Using Eq. (16), it is possible to reconstruct standard observables such as the forward-backward 



asymmetry, Ap-g, and the longtitudinal polarization fraction, Fi^: 

Afb = 1(231 + 3'^) and Fl = -S^. (19) 

o 

As explained in Sec. [T| our focus is on those observables that will be measurable at LHC6 
without a full-angular analysis. In order to keep the experimental complexity to a minimum, 
these observables should require information on only one or two of the angles. ApB; which 
depends only on 9i, and Plj which depends only on 9k*, are well known examples. They can 
be expressed as: 

d(r + r) /d(r + r)_ \_ 




/2 dq^ d cos 9k* / dq"^ 

where the latter expression makes use of the massless lepton approximation. We also study the 
possibility of an early measurement of 5*5, which can be measured using only cos 9k* and (p. It 
is possible to express this as 

dcoG^K- ^'(^"^^ / d(r + f) ^22) 
dq^ d cos 9k* dcp / dq^ 

A comprehensive study of the effects of the Wilson coefficients on the above observables, and 
vice-versa, can be found in Ref. [15]. We note that 3^, A^, and Ag can also be extracted 
by the counting of signal events over one or two angles. 3^ is related to the well known and 
theoretically clean observable Aj. [12]; to be precise, 3s equals ~ ™ massless 

lepton limit. While significant enhancement of a'^^ is possible in the presence of non-SM Cy*'^ 




|51j . the 2(1 — Pl) prefactor implies that the enhancement is less pronounced in 3^ [20]. The 
smallness of ^3 means that the experimental sensitivity to ('S'3)i_6Gev2 ^^^^ limited in the 
first few years of LHC6 data taking; thus, the study of 3^ is thus left for other works 
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Enhancements to Aj and Ag in the presence of NP phases can, however, be sizable [H] and 
could, in principle, lead to reasonable experimental resolutions, particularly for (^9)i_6GeV^' 
However, these measurements will still be experimentally challenging in the first few years. For 
these reasons we choose to focus on ^fb, -^l and for early study at LHC6. 

As stated earlier, NP enters the calculations through contributions to the Wilson coefficients; 
constraints on these contributions are described in the Sec. [H It is well known that for certain 
values of q'^, the observables ^fb and 6*5 vanish. We refer to these values of as the zero- 
crossing points, go(^FB) and (79(55). They are particularly sensitive to NP, and can be used 
to further constrain the values of the Wilson Coefficients. At leading order, in the large recoil 
limit, and for real values of the Wilson coefficients, it is possible to obtain simple expressions 
for [111 [31] andg2(55): 

go(^FB) = -2mBmb-^; q^iS^) = gcff ^ ^^(geff ^ g/eff) • (23) 

In deriving these results we make use of the soft form factors, following Refs |3H I41j. The two 
observables provide complementary sensitivity to NP through their differing dependence on the 
Wilson coefficients, and allow for sensitivity to both chiralities of O7. The cancellation of the 
soft form factors and the relative smallness of 0{as) corrections mean that both zero-crossing 
points meet the criteria for theoretical cleanliness given in, e.g., Ref. |13j . In addition, we define 
the gradient of ApB and ^5 at their zero-crossing points, 



dO 

dq' 



(24) 



where O is the observable ApB or S5 respectively, ^fb has also been studied in the context of 
B — )• KttI^I^ [52j, where expressions for qQ{ApB) and Go(j4fb) were determined for the case of 
an energetic kaon and soft pion. However, the kinematic region where the Ktt pair is energetic 
is dominated by the K* , and non- resonant effects can be neglected. 



3.2 Overview of Specific Models and Eff'ects on Wilson Coefficients 

The observables for — t- K*^fi^fx~ are most sensitive to the Wilson coefficients C!^^, Cg^, 
C^Q and their primed equivalents, so we concentrate on the NP contributions to these in this 
section. We also consider C^^ and Cp^ for completeness; however, experimental sensitivity to 
their effects are expected to be limited in this decay. 

• Flavour Blind MSSM (FBMSSM): Here the MFV version of the Minimal Supersym- 
metric Standard Model (MSSM) is modified by some flavour-conserving but CP-violating 
phases in the soft supersymmetry (SUSY) breaking trilinear couplings ^3j. The Wilson 
coefficients we use correspond to those calculated in scenario FBMSSM II defined in Table 
11 of Ref. [15]. The additional CP- violation contributes substantial complex phases to C|^, 
however there is no flavour structure beyond the SM, so primed operators are suppressed 
as in the SM. As in all SUSY models, scalar and pseudoscalar operators arise due to the 
additional Higgs doublet. 
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Figure 1: Theoretical predictions for vIfB) Fl and S^. The red (continuous) hne is the SM, the 
blue (dashed) line is the GMSSM, and the purple (dotted) hne is the FBMSSM. 



• General MSSM (GMSSM): Minimal flavour violation is not imposed, and generic flavour- 
and CP-violating soft SUSY- breaking terms are allowed [54J. The Wilson coefficients we 
use are close to the scenario GMSSM IV in Ref. [15], corresponding to large NP contri- 
butions to both Cj^ and Cy'^^ allowed by existing experimental bounds (see Sec. |4]). 

The Wilson coefficients in the above scenarios are given explicitly in Tab. [6| The central 
values for the distributions of A-pB-, Fh^ and are shown in Fig. [l] for the SM, the GMSSM, and 
FBMSSM, along with estimates of the theoretical uncertainties. The agreement with previous 
results is good. The predominant sources of the uncertainties are the form factors, hadronic 
parameters, and quark masses, which are determined as discussed in Sec. [2] We also include the 
uncertainty arising from varying the factorization scale, /i, in the range fj, G [/^/2, 2/i]. The three 
distributions all show significant variation for the models considered here, as do the position or 
absence of the zero-crossing points in Ayb and 5*5 in the range G [1, 6] GeV^. 

4 Constraints 

Experimental results can be used to constrain the NP contributions, denoted Cf^, to the Wilson 
coefficients: we define Ci = Cf^ + Cf^. We can then determine possible model-independent 
effects of NP on — )■ K*^ . The most important constraints on the Wilson coefficients 
are from the following measurements: 

• Branching Ratio for Bg — )■ fj-^ n^: This is used to constrain the possible NP contribution 
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Model 





SM 


FBMSSM 


GMSSM 


Cfifi) 


-0.306 


0.031+0.475i 


-0.186+0.002i 




-0.007 


0.008+0.003i 


0. 155+0. 160i 


Cf{ii) 


-0.159 


-0.085+0. 149i 


-0.062+0.004i 




-0.004 


-0.000+O.OOli 


0.330+0.336i 


ACf (/.) 


4.220 


4.257+O.OOOi 


4.231+O.OOOi 




0.000 


0.002+0.000i 


0.018+O.OOOi 




-4.093 


-4.063+O.OOOi 


-4.241+O.OOOi 




0.000 


0.004+O.OOOi 


0.003+0.003i 


(7s(/i)/GeV-i 


0.000 


-0.044-0.056i 


O.OOO+O.OOli 


Cp{fi)/Ge\-^ 


0.000 


0.043+0.054i 


0.001+O.OOli 



Table 6: NP Wilson coefficients at = mb,ps(2 GeV/c^) = 4.52GeV/c2 in the FBMSSM and 

"li where Cx{^^) = {Cx - C^)(^) iov X = S or P. 



GMSSM as described in Sec. 



to the scalar and pseudoscalar operators. To calculate the branching ratio we use the 
standard result from Ref. ITS] 



B{B., ^ = rBJlr^s.'%^\VM\^ l - 5^(|5p ( 1 - 5^ ] + |P|^), (25) 



mo V mi? 



with the definitions 



^ 2 

(Cs - C's); P='^{Cp- C'p) + m,[Cl^ - C[f). (26) 



We use /s, = 0.259 ± 0.032 GeV [55j, tb, = 1.456 ± 0.03ps [56] and m^^ = 5.37 GeV 
|43j . and other numerical parameters as in Ref. [15j . In agreement with existing results, 
we find the SM prediction, BR{Bs pi'^l^') = (3.70 ± 0.31) • 10"^ to be well below the 
current experimental upper bound 3.6 • 10~® |57| . 

Branching Ratio for B — Xgl^l^: We compare NP predictions for B{B — Xsl^l^)i_QQgy2 
to the mean experimental value (1.60 + 0.51) • 10~^, as adopted in Ref. jT3|, combining the 
results of BAR4R, (1.8 + 0.7 + 0.5) • 10"^ 158], and Belle, (1.49l[!:^^ + 0.50)-10-6 [59]. This 
helps to constrain the NP contribution to Cjg^iQ as well as C^p. As an inclusive mode, 
the calculation for the region G [1, 6] GeV^ of the branching ratio is theoretically clean. 
We use the expression for the differential decay distribution in Ref. [60], but also include 
the NLO corrections computed in Ref. [61] , and the contribution of the primed operators 
as in Ref. [62]. Using our parameters we find B{B — )• Xsl^l^) = (1.96 + 0.11) • 10~^ for 
the SM. 

Branching Ratio for B — )■ Xs'j- The current experimental average for > 1.6 GeV is 
B{B Xs-f) = (3.52 + 0.23 + 0.09) • 10"^, as calculated by the Heavy Flavor Averaging 
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Observable 


Wilson Coefficients 


^FB 




rics 


cf 


Fl 


r eff /-^/ off 




S5 


rcff 
U7 , 


r/eff 
<-^7 ' 


reff /^/ efF 
, Oj^o 



Table 7: Relevant observables and the Wilson coefficients they most strongly depend on |15| . 

Group [56J. We use the recent theoretical SM result of Ref. [63], (3.28 ± 0.25) • 10"^ for 
> 1.6 GeV, and include NP effects as in Ref. [51]. The SM calculation makes use of 
the kinetic renormalization scheme for determining rric and nib; an alternative calculation 
using the IS scheme leads to a branching ratio of (3.15 it 0.23) • 10~^ |64| I65j : however, 
our results are not sensitive to the difference between these two values. 

• Time dependent CP Asymmetry S{B — )• K*^): This constraint is sensitive to the 
photon polarization, and, hence, to Cj^^ . This should be compared to S{B — K*^) = 
(-1.6 ± 2.2) • 10-1 from experiment [56]. Our SM result S{B K*-y) = (-0.26 ± 
0.05) • 10~i agrees with that of Ref. [H] within uncertainties. In Refs [66l EZj, the soft 
gluon contribution was calculated, leading to a small correction to the predicted value. 
This is neglected in our treatment as it has little effect on the constraining power of the 
experimental measurement. 

• Integrated Forward-Backward Asymmetry (^fb)i.6GcV^ — ^ K*^fi^fi^ : 
We use the existing — )• K*^ ji^ jjL~ measurements as constraints. Recently Belle has 
made a measurement of the forward-backward asymmetry, and finds the integrated A-p-Q 
value in the region 1-6 GeV^ to be -0.26 ± 0.29 [8]. This is to be compared to our SM 
prediction of 0.04 it 0.03, which is in agreement with the recent result in Ref. [68] • This 
observable constrains the Wilson coefficients as seen in Tab. [71 We look forward to a 
1-6 GeV^ measurement from CDF with great interest [S]. 

• Integrated Longitudinal Polarization Fraction (^L)i-6Gev2 ^d, K*^ ji^ : 
Belle has also recently measured the Longitudinal Polarization Fraction to be 0.67 it 0.24 
[S]. This should be compared to our SM prediction 0.76 it 0.08, also in agreement with 
Ref. [HS]. Again this constraint affects Wilson coefficients as seen in Tab. [7| 

In order to assess the impact of these constraints on the NP contributions to the underlying 
Wilson coefficients in as general a way as possible, we have performed a semi-random walk 
through parameter space. We allow {Cs — Cg), {Cp — C'p) and the NP components of Cy'lio 
vary simultaneously, both in magnitude and phase. To our knowledge this has not been done in 
previous studies. At each randomly chosen point in parameter space, predictions are made for 
the six observables listed above. The point is then either accepted or rejected using a modified 

metric that treats experimental uncertainties as being normally distributed, but theoretical 
uncertainties as having uniform probability within the specified range. Following traditional 
minimization techniques, the random walk is guided by this modified so that regions with 
lower values may be identified. Using this method, a sample of 2.5 • 10^ independent sets of 
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Observable 


Experiment 


SM Theory 




3.6-10-8 [57] 


(3.70 ±0.31) - 10-9 


B{B ^ XsI+I-UgcY^ 


(1.60 ±0.51) - 10-6 


(1.97 ±0.11) - 10-6 


B{B ^ X,7) 


(3.52 ±0.23 ±0.09) - 10-^ [56J 


(3.28 ±0.25) - 10-^ 


S{B K*j) 


(-1.6 ±2.2) - 10-^ [56] 


(-0.26 ±0.05) - 10-1 


(^FB)i_6GeV2 


-0.26 ±0.29 [H] 


0.04 ±0.03 


(^l)i-6GcV2 


0.67 ±0.24 [H] 


0.76 ± 0.08 



Table 8: Experimental measurements used as constraints, along with theoretical predictions in 
the SM. 



Wilson coefficients was produced. Each set results in predictions for the observables listed above 
with better than 2a agreement with current measurements. It was found that the agreement 
between existing measurements and the SM is excellent, with a P^r degree of freedom of 0.35. 
While this is not implausible for six degrees of freedom, the level of agreement suggests that 
more detailed study of the theoretical uncertainties will be required as experimental resolutions 
improve. 

Fig. [2] shows the range of values found for the phase and magnitude of the NP contribution to 
Cy^ and Cj^^ (at the scale fj, = mw) during the parameter space exploration. The colour index 
shows the mean value of the probability that a point is compatible with current experimental 
results. Areas with probability greater than la are shaded red, while those with less than la 
are shaded blue. The outline of the la contour can clearly be seen. The values of the Wilson 
coefficients for the SM, FBMSSM, and GMSSM are also shown. 

Fig. [2] can be compared to Fig. 2 from Ref. [H], in which and C^'^^ are assumed to be 
real and all other Wilson coefficients SM-like. The effects of weakening these assumptions can 
be seen. Similar figures are shown for the other Wilson coefficients in Figs [3] and |4j The allowed 
regions of parameter space are still large, particularly if NP phases are allowed. In contrast to 
Ref. [14J, constraints from measurements at high-g-^ (low recoil) are not included as we 
feel that NLO effects are not under control in this region. The effect of this constraint may be 
seen by comparing our Cfg figure, shown in Fig. [Zj with that in Fig. 2 of Ref. [H] . 

The ensemble of constrained NP models can also be used to explore the likely values of the 
A-pB and ^5 zero-crossing points in the range q"^ G [0.5, 15] GeV^. While it should be noted that 
theoretical uncertainties are not well controlled over this range, the majority of points within 
the la contour lie within the theoretically clean region, G [l,6]GeV^ (see Fig. 5a). It was 
found that 8% of the parameter space points considered had no ApB zero-crossing in the range 
q'^ G [0.5, 15] GeV^. For S5, only 2% of points had no zero-crossing in the same range. Fig 



5b 



shows the Ape and ^5 gradients at their zero-crossing points. We find that, for the majority of 
points, ^0(^5) is greater than Go(Afb)- This will have an impact for the ^gC'^s) experimental 
analysis discussed in the Sec. 5.4.2[ 

To summarize, in this section we have considered six existing experimental constraints, and 
used these to determine the allowed regions in parameter space for the NP contribution to 
the Wilson coefficients. These allowed values for the Wilson coefficients were then used to 
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Figure 2: [Colour online] Allowed parameter space for the NP contribution to and C^^^ 
at the scale = mw, as described in Sec. |4} Points with a compatibility with data of 68% or 
better are drawn with a dark (red) colour palette, while those with less than this are drawn 
with a light (blue) palette. The SM point is shown in black at the origin, while the FBMSSM 
is a green square and the GMSSM is a blue triangle. The Wilson coefficients for these models 
are shown in Tab. |6j 

find corresponding predictions for ggl'^s)' 9o(^fb)5 GoiS^), and Go(^fb)- In the following 
sections, we investigate the experimental sensitivity to the observables ApB, S^, and Fl, and 
how measurements of these could have an impact on the allowed NP contributions to the Wilson 
coefficients. 



5 Experimental Sensitivities 



Three observables that can be measured as a function of by counting signal events in specific 



angular bins, using Eqs (20)-(22), were highlighted in Sec. 3.1: ^fBj -^L) and ^5. These 
observables should be suitable for early measurement at LHC6. In the following, we estimate 
the experimental sensitivities in order to make a fair comparison between these observables. 

LHC6 is expected to collect ~ 6.2 • 10^ signal events per 2fb~^ of integrated luminosity 
with a signal to background ratio of approximately four [69^ I70j . With relatively small data 
sets it should be possible to extract the values of these observables integrated over q^. These 
measurements provide an early opportunity to discover NP in 6 — t- s transitions. For larger 
data sets it will be possible to map out the dependence on q'^ as well, allowing for additional NP 
discrimination. Studies of these two approaches can be found in Refs [181 ttl] for tlie observable 

To assess the impact of each potential measurement on the allowed NP parameter space, 
simple analyses have been developed to extract the integrated values of ^fb, -^l, and ^5 
in the region q"^ E [l,6]GeV^. In addition, analyses have been constructed to extract the 
q^ dependence of ^fb and ^5, along with their zero-crossing points; the latter can be found 
numerically from the Apb{q^) and /S'5(g^) distributions. In order to minimize the experimental 



16 



-3-2-10 12 3 -3-2-10 12 3 




Figure 3: [Colour online] Allowed parameter space for the Wilson coefficients C^lx^ after ap- 
plying relevant b ^ s experimental constraints. The colour coding is the same as in Fig. [21 
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arg[Cs-Cs') argtCp-Cp) 

Figure 4: [Colour online] Allowed parameter space for the Wilson coefficients c[q'^^ and {Cs,p — 
C'g p) after applying relevant 6 — >• s experimental constraints. The colour coding is the same as 
in Fig. [2] 
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10 12 14 



(a) 




Figure 5: [Colour online] Fig. (a) shows allowed values of the ^fb and ^5 zero-crossing points 
in the range G [0.5, 15] GeV . The SM point and its uncertainty is shown as a black ellipse. 
Fig. |(b)| shows the gradient of the ^fb and at the zero-point. For comparison, the line 
^0(5*5) = Go{Afb) is included. In each case the colour index has the same meaning as in Fig.[2j 



uncertainties on these points, a larger region of G [0.5, 8.5] GeV^ was used for these analyses 
following Ref. [19]. An ensemble of 1200 simulated Bd — K*^{-^ K7r)fi'^fj.~ data sets was 
created, each containing the (Poisson fluctuated) number of signal and background expected 
from 2fb~^ of integrated luminosity at LHC6. Other integrated luminosities were obtained by 
linearly scaling the yield estimates. Each analysis was then run in turn on the data sets in order 
to estimate the statistical uncertainty expected for each measurement. This allows for a fair 
comparison to be made between observables for a given integrated luminosity. 



5.1 Data Set Generation 

The theoretical framework introduced in Sec. [2] was implemented as a plug- in for the standard 
decay tree simulation tool EvtGen |23J. This allows Bd — > K*^{^ K7r)//+/i~ events to be 
simulated. A simplified background sample was generated separately. This was flat in the 



three decay angles defined in Sec. 3.1 but followed the signal distribution in q and a gently 



falling exponential in the B invariant mass, uib- All events had tub within a wide window 
^^5Q MeV around the nominal B mass. A central signal region was also defined with width 
±50 MeV. Events outside of this region were assumed to be part of a background dominated 
side-band. Signal and background events were generated following the relative normalization 
given in Refs [69l [TO]. For each event in a data set, the three decay angles, q'^ and tub were 
determined and used as input for each analysis. 



5.2 Integrated Analyses 

The integrated quantities can be extracted by estimating the number of signal events in each 
angular bin using a fit to the tub distribution. The signal contribution was parametrized as a 
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Gaussian with an exponential tail, while the background was modelled as an exponential with a 
negative coefficient. A fit was performed to each data set to extract the signal and background 
shape parameters for that sample. Each sample was then reduced into the relevant angular bins. 
For ^FB) following Eq. (20) these bins would be cos^; G cos 6i £ [0, 1] for all events 

in the range G [1)6] GeV^. To extract an estimate of the number of signal and background 
events in each angular bin, a separate fit to the itib signal and background distributions was 
then performed, keeping all shape parameters fixed. The value of (^FB)i_6Gev2 determined 
with Eq. (20). A similar procedure was applied to Eqs (21) and (22) to extract (-^L)i_6Gev2 

('S'5)i_6GeV2- 



5.3 Dependent Analyses 

Following Ref. [19], a polynomial shape was fit to the distribution in each angular bin. The 
method proceeds as in Sec. |5.2[ using the B mass distribution to find the total number of signal 
and background events in each angular bin. However, the background shape extracted is used 
to estimate the number of signal events in the B mass signal window. The dependence of the 
signal and background distributions was parametrized using second and third order Chebyshev 
polynomials respectively. A simultaneous fit in the signal and side-band regions of the B mass 
distribution was used to determine the shape parameters of signal and background polynomials 
using the relative signal/background normalization found from the B mass fits. In the case 
of ^FB, the procedure would lead to the extraction of two q'^ dependent signal polynomials: 
one for events with cos0; G [~1)0] the other for cosOi G [0,1]. The value of ApB {q'^) 



can then be found using these polynomials and Eq. (20). The ApB zero-crossing point was 
found numerically from the combined functions. A similar approach was applied to and its 
zero-crossing; however, six angular bins in 9k* and (j) were required. 



5.4 Results 

When comparing different observables and analyses it is useful to consider the mean expected 
experimental sensitivity for a given integrated luminosity. These expected sensitivities can be 
calculated from the ensemble of toy LHC6 experiments introduced in Sees |5.2| and |5.3] 1200 
individual experiments were performed, and for each one a value of, for example, gQ(^FB) was 
found. Following Ref. [13], the mean, one and two sigma contours could then be found from these 
results. The method used allows for non-normally distributed results by putting the ensemble 
in numerical order and then selecting the values closest to the contoui|^ Any biases introduced 
can be identified by comparing the median result and input value. Example ensembles are 
shown in Fig. [t] for qQ{ApB) and (7q(55), assuming 2fb-i of LHC6 data and following the SM. 

^For the one sigma bound these would be the 188*'^ and 1010**^ resuhs in the ordered ensemble. 
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2 fb~^ 


1 fb^^ 


n 5 fb~^ 


1 A \ 

(^FB)i_6GcV2 


+0.03 
-0.04 


+0.05 
-0.03 


+0.08 
-0.06 


/ 77V \ 


+0.02 
-0.02 


+0.04 
-0.03 


+0.04 
-0.06 


('S'5)i_6 GeV^ 


+0.07 
-0.08 


+0.09 
-0.11 


+0.16 
-0.15 


9o(^fb) 


+0.56 
-0.94 


+1.27 
-0.97 






+0.27 
-0.25 


+0.53 
-0.40 





Table 9: Estimated 1<t LHC 6 sensitivities for 2 fb ,lfb and 0.5 fb of integrated luminosity, 
assuming the SM. 




2 4 6 8 2 4 6 8 



(GeV^) (GeV^) 

Figure 6: Projected experimental sensitivities to the observables ^fb and 5*5 using an unbinned 
polynomial fit to 2 fb~^ of LHC6 data in the range q"^ G [0.5, 8.5] GeV^. The dashed line shows 
the input distribution, while the solid line shows the median of an ensemble of 1200 fits. The 
light and dark contours show the estimated one and two a contours. 

5.4.1 Integrated Quantities 

The estimated la sensitivities for the integrated observables (^FB)i_6Gev2) (^L)i_6Gev2 ^^^"^ 
(5'5)i_6Gev2 LHC6 data set sizes of 2fb~^, 1 fb~^ and 0.5 fb~^ are shown in Tab. [9| Any 

differences between the input and extracted median values were seen to be small relative to 
the estimated uncertainties. The estimated LHC6 experimental uncertainties are of a similar 
size to the current theoretical uncertainties, and much smaller than the current experimental 
constraints [8j. 

5.4.2 Zero-Crossings 

Fig. |6] shows the projected experimental sensitivity to the full ApB and 5*5 distributions for 
2fh-^ of LHC6 SM data. For ease of comparison with SM predictions, the zero-crossing point 
is extracted from the dependent distributions. These are shown in Fig. [t] for the same data 
sets as used in Fig. [6j The estimated la uncertainties are shown in Tab. [9} As discussed in 
Ref. |17] . the experimental uncertainty will scale approximately linearly with the gradient at 
the zero-crossing, leading to the large difference in estimated sensitivities seen for qQ{ApB) and 
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Figure 7: Projected experimental sensitivities to the zero-crossings of A-p-Q and 5*5 using an 
unbinned polynomial fit to 2fb~^ of LHC6 data in the range G [0.5,8.5] GeV^. The colour 
coding is the same as in Fig. l6j 



ql{S^) inTab.[9| 

The difference in gradients between ^fb and ^5, seen in Fig. 5b for the majority of NP 
points, makes q^iS^) an attractive experimental target, assuming that any practical difficulties 
associated with the 6k* and decay angles can be overcome. We see that the relative steepness 
of the 55 distribution is such that the experimental uncertainty on qQ{S^) should be compet- 
itive with that on (j^^Ayb) for the majority of the allowed regions of parameter space. For 
0.5 fb^^, biases on the zero-crossing points become significant when using the unbinned analysis 
technique; however, it is likely that coarse estimates of gQ(ylFB) and (70(5*5) could be extracted 
even at this relatively small integrated luminosity using alternative techniques, such as those 
discussed in Ref. [181 ■ 



6 Impact of Future Measurements 

The relative impact of the different analyses presented in Sec. [5] can be assessed by revisiting 
the parameter space exploration performed in Sec. |4j We are interested in how including these 
new measurements would affect the current constraints on parameter space. It is assumed that 
LHC& will make 2fb ^ measurements of the observables (^FB)i_6Gev2> {S^)i-&G(N^i 1-6 GeV^ ' 
?o(^fb)i and q^iS^) and that the resulting experimental uncertainties are symmetrized versions 
of those given in Tab.[9j In addition, we assume that the measured values of these observables are 
not affected by NP, and are as given in Tab. [sj The total for each point in parameter space is 
then updated to reflect these hypothetical SM measurements. Where individual measurements 
are superseded by LIIC6 measurements, they are replaced with no attempt at combination. 
However, other constraints, such as B{B — )• X^^), are included as before. In this way the 
constraining power of each analysis can be compared. 

Fig. [s] shows the relative impact of these measurements on the NP component of C^^. In 
Fig. ISal SM values of (^FB)i_6Gev2 ^^"^ 9o(^fb) are imposed with the estimated 2fb~^ experi- 



mental sensitivities taken from Tab. M Fig. 8b shows the impact of (-FL)i_6Gev2' '^liile Fig. 



8c 



shows the impact of both (5'5)^_gQgY2 and q^iS^) for the same LHC6 integrated luminosity. 
These should be compared with the currently allowed Cj^ parameter space shown in Fig. [2j 
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Figure 8: [Colour online] The relative impact of different proposed LHC6 measurements after 
2fb~^ of integrated luminosity, assuming the SM, on the NP component of Cy^. In each case 
the colour index has the same meaning as in Fig. [21 
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Figure 9: [Colour online] Allowed parameter space for the Wilson coefficients Cf, C'^"^, Cf 
and C[f^ after 2fb~-^ measurements at LHC6 of {Fi^) i_q Qey2 , (AFB)i_6Gev2' 9o(^fb), (5'5)i_6Gev2 
and 00(85), assuming the SM. The colour coding is the same as in Fig. pi 



The small statistical uncertainty found in Sec. [5] for (/qC'S's) provides a stringent constraint on 
parameter space. This emphasizes the importance of an early measurement of S5, in addition 
to ^FB and Fl. 

Fig. [9] shows the combined effect of the measurement of the proposed observables, again 
assuming the SM and the estimated sensitivities from Tab. [9] for the NP contribution to the 
Wilson coefficients Cf^, Cj'^^, Cg^ and C[q^. The amount of parameter space left after these 
measurements would be significantly reduced, with most NP contributions excluded at the la 
level unless there are large NP phases present. This again illustrates the importance of CP 



observables as described in |14| I15j. The FBMSSM and GMSSM models from Sec. 3.2 could 
also be excluded at better than 95% confidence in this case. 
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7 Summary 



A new next-to-leading order model of the decay -B^ — )• K*^n^fi^, that features QCD factoriza- 
tion corrections and full LCSR form factors, was presented. This includes an expression for the 
decay amplitude in terms of an updated set of auxiliary functions; these can be compared dir- 
ectly to the previous model, based on Ref. |22j . The auxiliary functions have been extended to 
include the effects of primed, scalar, and pseudoscalar operators, which may become important 
in certain NP scenarios. 

The observables ^fBj -^Li and S5 were identified as being promising for a relatively early 
measurement at the LHC, as they can be extracted as a function of by counting signal 
events in specific angular bins, using Eqs (20)-(22), and correspond to large features in the 
angular distribution. We also obtained a simple expression for ggC'S's) at leading order, in terms 
of Cj^, C'y'^^, and Cg^, and showed that it has reduced hadronic form factor uncertainties in 
the large- recoil limit. Considering current experimental constraints leads to restrictions on the 
possible NP contributions to the Wilson coefficients. The allowed values of the ^fb and 5*5 
zero-crossing points, and the gradient of the ApB and 55 distributions at these points, were 
explored. The relative steepness of the distribution, even in the presence of NP, makes 
q^iS^) an experimentally attractive target, as it will lead to a smaller experimental uncertainty. 

In order to investigate the impact of measuring the proposed observables on the NP con- 
tributions to the Wilson coefficients, and to compare their relative impact, we estimated their 
sensitivities at LHC 6. We studied the sensitivity to the q^ integrated values and zero-crossing 
points of ApBi Fl, and 5*5. The prospect of measuring S5 and its zero-crossing at LHC 6 has 
not been previously explored. 

Using a combination of (i^L>i_6Gev2' (^FB>i_6Gev2' 5o(^fb), (5'5>i_6Gev2' and q^iS^), we 
showed that 2fb~^ of LHC6 data could greatly reduce the allowed parameter space. The 
contribution of ^5 to this is very significant and can, in part, be attributed to the small statistical 
uncertainty expected on (7g(S'5). We have also shown that if the decay is SM-like, the GMSSM 
and FBMSSM points considered would be ruled out by LHC6 with 2fb^^. We conclude by 
stressing that making measurements of 5*5 and its zero-crossing would provide an interesting 
and complementary measurement to others currently planned. Ba — )■ K*^fj.'^^~ is a promising 
channel for constraining models or making a NP discovery. We look forward to the ffist LHC 
results for this decay. 
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Appendix 



A Operator Basis 

The effective Hamiltonian for — )• K*^fi^fi~ can be expressed in terms of effective operators 



and Wilson coefficients as described in Sec. 2.2 We provide explicit expressions for a subset of 
these operators, which play a key role in the decay. Definitions for the remaining operators can 
be found in Ref. 



(27) 
(28) 
(29) 
(30) 
(31) 
(32) 



where g is the strong coupling constant, e is the electron charge, is the b quark mass in the 
MS scheme, as described in Sec. 2l and Pl,r = (1 =F 75)/2- 



Or = 


gZ 


0'r = 


^mi,{sa^,PLh)F^'\ 

gZ 


Os = 


9 




-7fife(sa^,r'^Pi&)G^-« 

9 


o,= 


e2 

gZ 


o',= 


e2 

^(s7/.^i?&)(A7^/^), 
9 


Oio = 


gZ 




e2 
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^Q^2^b{sPRh){^l^l), 


0's = 




Op = 


mb{sPRb){fi-f5n), 


o'p = 





B Angular Coefficients 



3.1 



and 



Here we provide the relations between the angular coefficients, l\ \ defined in Sec. 

I~l (sic) 

the auxiliary functions defined in Eq. (8). We first express the II s in terms of transversity 
amplitudes as in Ref. [jl 5J . 
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4m, 



+ 



At\-' + 2ReiA(^A§^) + l^MsW 
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|Ai|2 + |Af|2 + (L^i?) 
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I A 
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\A^\^ + iL^R) 



V2 
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(38) 
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-'6 
-'6 

h 
Is 



2/3 
4/3- 



m 



m,. 



Re [A^A*s + (L ^ i?)] 



\/2/3 
1 



lm{A^A 



m 



/3^ Im(.4^Ai*) + (L ^ i?) 



/g = rim(ArAi) + {L^R) 



(39) 

(40) 
(41) 

(42) 

(43) 
(44) 



These transversity amplitudes are projections of the decay amphtude onto various combinations 
of hehcity states of the K* and the virtual gauge boson. The projections can be achieved by 

1/2 

contracting 7^ with the virtual gauge boson polarization vector. We use four basis vectors for 
the virtual gauge boson polarization vector corresponding to transverse (it), longtitudinal (0) 
and time-like (t) states, and three basis vectors for the virtual gauge boson polarization vector 
corresponding to transverse (it) and longtitudinal (0) states. One first extracts the helicity 
amplitudes and Hq using the basis polarization vectors +,-)0 respectively for both the 

K* and the virtual gauge boson. Ht is found by taking the longtitudinal polarization vector for 
the K* and the time-like polarization vector for the virtual gauge boson. Using the relations 



■/II 



V2 



(45) 



and ^0 = Ho, At = Ht, one then obtains expressions for the transversity amplitudes in terms 
of ^(g2) to 5(g2), 



N ruB 



1 — m 



K* 



ai{q^) + Xh{q^ 



niK* \/q- 

"^^^^ {F{q') - (1 - m^.)G(g^) - eH{q')) 



(46) 
(47) 

(48) 
(49) 



TUK* V9 

where i = L/R. We use the standard normalization and definitions following Ref. |12j . 



'1 



4m?, 



A =1 + m^, + - 2 {f + ml., (l + f)) 



(50) 
(51) 
(52) 
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where a is the electromagnetic couphng constant and Gi? is the Fermi constant. In the above 
definitions of the transversity amplitudes, the functions ^L/i?^^^), ci/ji{q'^), are ana- 

logous to those defined in Ref. [7T] . 



aL/R{q^)=B{q^)TF{q^), 



(53) 



(54) 



(55) 



Using the above it is possible to compare the predictions of Eqs ([8]) to the standard results in 
the literature, and we agree with Ref. |15j . 
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